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Abstract. We have carried out a thorough optical study of the closest star formation jet. The inner arcminute surrounding the 
class O/I binary protostar L1551 IRS5 and its associated jet (HH154) have been observed using the Hubble Space Telescope 
with the WPPC2 camera, the ESO New Technology Telescope with the EMMI spectrograph and the Nordic Optical Telescope 
with the ALFOSC spectrograph. This data set is compared to earlier ground based imaging with the aim to study the evolution 
of this particular jet, and its possible interaction with the molecular material in the bipolar molecular outflow associated with 
this source. The velocity field of the jet is mapped out. The highest velocities are found in the vicinity of the recently discovered 
X-ray source emanating from a shock in this jet. The energy radiated by the X-ray source is compatible with these velocities. 
The Ha and H/3 emission from the jet is used to determine the extinction, which is found to increase inwards in the jet towards 
the protostar. The extinction towards the X-ray source is consistent with the one determined from the X-ray spectrum. 
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pre-main sequence 


1. Introduction 

Among bipolar molecular outflows, the one centered on 


the youn 
covered 

g stellar object IRS5 in L1551 was the first dis- 

Sne 

1 et alJ 

1980t). It is relativelv nearbv, ~150pc 

jKenvon et al. 

119941 

Reiourthlll999ll. allowing the studv of 


small spatial elements, and oriented with its major outflow axis 
most li kely at an angle of between 45° and 60° dUiseau et alJ 
l2005bl) with respect to the plane of the sky. This results in 
red- and blue-shifted CO outflow lobes which are spatially well 
separated, and which have an angular extent of several tens 
of arcminutes on the sky. The red- and the blue-shifted CO 
lobes contains several Herbig Haro objects (HH objects) most 
pr ominently seen in the blue- shift ed lobe, (e.g. HH28, HH29). 
In ICudworth & Herbid ( ll979l) and iFridlund et all ill 984 fl^ 
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* Based on observations made with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica de Canarias, and 
based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. 


Il998l> the authors argue, based on proper motion vectors, ve¬ 
locity gradients in CO, and the radial velocity held, that the 
outflow is originating at IRS5 and is interacting with the am- 
bient medium at the pos tion of the Herbig Haro objects. In 
iDevine et alJlll999ll2000l) . the authors on the other hand argue 
that the HH28 and HH29 are excited by a flow from another 
YSO, i.e. L1551NE. 

The low mass characteristi cs of IRS5 were dete rmined by 
observations in the far infrared (iFridlund et aljl9 ^. A jet has 
been observed at optical waveleng ths (iMundt & Fried ll983[ 

iNeckel & Staudelll987t ^idlund & Liseaulll994nT99^ in the 

near IR dcampbell et^l. 1988l:IPvo et alJ2002il2005l) an d in th e 
radio continuum dCohen et al .I ll982t iRodriguez et al.lll984 . 
The realisation that th is jet has Herbig Haro characteristics 
(ISarcander et aljfl985ll has le d to it receiving the designation 
HH154 in the compilation of iReipurthI (Il999ll. Morphologica l 
cha nges were detected within the jet bv ISarcander et al.l lll985l) 
and iNeckel & Staudel (Il987l) . The latter found that over the 
epoch 1983 - 1987, a new knot appeared near IRS5, and that 
the transverse velocities for several knots in the jet were of or¬ 
der ~ 190 km s ' - in agreement wit h earlier proper motio n 
measurements of one of the knots by ISarcander et alJ lll985l) . 
INeckel & Staudel (Il987l) introduced a consistent nomencla¬ 
ture for the features observed along the jet. This was also 
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used by iFridlund & Liseaul ll 19941) - hereafter FL94 and 
Fridlund & Liseau lll998l) - hereafter FL98, and we continue 


to use their designation in this paper, consistently updated in 
view of the evolution of the jet. 

At all observed wavelengths, the jet is relatively compact 
compared to the molecular outflow, growing from ~ 7 arc- 
sec to ~14 arcsec between 1983 and 2001. The jet terminates 
in a bright emission knot - ‘D’ according to the nomencla¬ 
ture described above - although some fainter, isolated emission 
peaks (A, B and C ) are also detected further ‘downstream’. 
iMundt et al.l lll99ll) used image restoration techniques on a 
deep [S It] image and suggested that the jet consists of a limb 
brightened cavity, with a high velocity wind permeating the 
cavity, which through interactions with the walls of the cavity, 
produce the HH knots. The optical properties of the jet were 
then further elucidated by FL94 who reported on high spatial 
resolution monitoring observations of the L1551 IRS5 jet car¬ 
ried out between 1989 and 1993, made from the ground in ex- 
cellent seeing conditions . FL94 confirmed the earlier results by 
Neckel & Staudej l[l987l) that morphological changes along the 


jet have dynamical time-scales of < a few years, and that these 
changes can be interpreted as proper motions of knots within 
the jet. Disappearance of ‘old’ and appearance of ‘new’ knots 
was recorded. It was also found that the brightest feature (D) 
in the jet accounts for more than 50% of the observed total 
emission in the lines of Flo' and [Sii] 6716,6731 A. Also from 
these observations, three different transversal velocity compo¬ 
nents (120 kms^', 190 kms"' and > 310 kms"') were identi¬ 
fied. Taken together with earlier radial velocity measurements 
found in the literature, FL94 hypothesized that the spatial dis¬ 
tribution of the knots was best represented by the features mov¬ 
ing close to the surface of a conical volume. FL98, which re¬ 
ports preliminary results from the first visit in the two-epoch 
Hubble Space Telescope (HST) program fully described in the 
current paper, interprets their data together with some spectro¬ 
scopic data as indicative of the existence of two separate jet- 
flows originating from IRS5. Each of these is possibly orig- 
inating from one o f the two components of IRS5 found by 
iLoonev et alJ lll997l) . FL98 also identified the brightest compo¬ 
nent (‘D’) of the jet with the ‘working surface’ of the ‘northern 
jet’ a nd discussed the structure of that shock. Using the mod¬ 
els of iHartiganl ill 9891) . they found that the ionisation fraction 
appears high (~ 50% - 100%). From the velocities determined, 
an estimate of the forces required to accelerate the individual 
knots could be made. This was compared to the force required 
to accelerate the molecular gas in the same volume of space as 
the jet (Fridlund & Knee, 1993) and found to be at least two 
orders of magnitude too low assuming an ionisation fraction of 
~ unity. FL98 pointed out the impact of this result on the ques¬ 
tion of the level of ionisation in optical jets which has turned 
into a major problem when trying to unify models for optically 
visible jets with those of the accompanying molecular outflows. 

The structures closest to IRS5 have been observed with 


interferometers in the radio contini 
GHz), by a number of researchers 

mm (1.4 GHz to -220 

iBieging & Cohen 1 19851 

Snell et alJ 198^ Rodriguez et alJ 

1986 

Keene & MassonI 

1990t Ladd et alil995l Looney et alJ 

mi 

). The latter authors 


pointed out that a 3 component model emerges with a large 


scale (10 to 30 arcsec) envelope, a disk like feature (~ 1 arc- 
sec) and a binary source with individual components separated 
by 0.35 arcsec, and where the disks have each a FWHM of 
^ 0.3 arcsec. This double source, which was originally dis¬ 


covered b\jBi£ging_^ Coheij 11^ and at that time inter- 


I 


preted by iRodriguez et ^ (Il986h as the inner ionized edges 
of a gas and dust t oroid, is now seen as a yo ung binary system 
(also suggest ed bviBieging & Cohen Il98 5l) with a separation 
of ~ 50AU llToonev et alJ T^tI) . llRodnguez et al ] l2003albl) 
have analysed interferometric radio continuum observations, 
and trace the two jets to their origin. Based on that data, those 
authors suggest that the northern jet-component emanates from 
the sou thern binary and vice versa. It has also been demon¬ 
strated jFridlund et al.lI2flfl2t Iwhite etlfl ]B) that there is 
also a fourth component in the form of a large (~ 15 000 AU), 
massive (several Mq), rotating flattened envelope of molecular 
gas. 


In the L1551 cloud, we And the astrophysical jets located 
nearest to the Earth (HH154, HH30 and the HE Tau jet). 
This implies that for HH154, the HST can achieve similar 
(~ 10'^^ cm), or better, spatial resolution t o that used in cur rent 
2D hydro-dynamical calculations of jets jRaga et alJl2004 . A 
multi-cycle observing programme with the HST consequently 
appeared well motivated, with the primary goal being to study 
the time evolution of the jet, to use the shock diagnostic emis¬ 
sion lines of Her and the [S ii] doublet at AAblll, 6731A in or¬ 
der to determine its physical properties. To detect fainter details 
we also obtained R-band (E675W Alter) observations. Eurther 
we wished to study the structure furthest in towards IRS5. 
Since the jet here suddenly appears from behind ~ 100 mag¬ 
nitudes of obscuring material (presumably the accretion disk), 
we try to observe light scattered up along the major axis of 
the jet. To achieve this we selected the E814W I-band Alter. 
As mentioned above, preliminary results have been reported in 
IFridlund et al.l (119971) and FL98. Included in the latter paper 
are some (also preliminary) ground based spectroscopic data, 
demonstrating a radial velocity gradient along the jet, as well 
as two separate velocity fields. 

Rece ntly, X-rays eman ating from the jet haye been dis- 
coyered (iFayata et alJ2002l) . The obseryations carried out with 
XMM-Newton haye a positional uncertainty of roughly the 
same size as the jet. The X-ray source with a temperature of 
a few million degrees and a yisual extinction of ~ 5 magni¬ 
tudes, was tentatiyely identified as originating from the work¬ 
ing surface (’knot D’), mainly because of its brightness, and 
feature D be ing the only source with reported lOiiil emi ssion 
then known dCohen & Fulleill98.4LSarcanderet al.ll 98.51) . The 
X-rays had to originate in the jet, since the star itself is hid¬ 
den by the yery large extinction from the surrounding disk. 
Th e position o f the X -ray source haye recently been measured 
bv iRallv et all ll20n3l) . who obseryed with the Chandra X-ray 
obseryatory proyiding a spatial resolution of ~ 1 arcsec. This 
has enabled them to determine that the X-rays emanate roughly 
from ~ 0.5 arcsec to 1 arcsec from the radio position of IRS5. 

In the present paper we giye the detailed results and inter¬ 
pretation from two HST epochs (cycle 5 & cycle 6). We also 
use new, extensiye ground based imaging and spectroscopic 
data to support our interpretation. 
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We confirm the dual structure of the jet itself, and we de¬ 
scribe the somewhat surprising morphology found close to the 
base of the jet. The time evolution of the jet is described, as dis¬ 
cerned from the measured proper motions of individual knots 
within the jet. We also discuss the structure of different indi¬ 
vidual shocks, including the excitation conditions and the elec¬ 
tron density along the jet and within separate knots, identify 
the velocity component associated with the X-ray source and 
the implications of having a strong X-ray source illuminating 
the accretion disk surrounding 1RS5. 

2. Observations and data reduction 

2.1. HST Imaging 

The area surrounding the LI551IRS5 jet was imaged using the 
Wide Field and Planetary Camera 2 (WFPC2 - Trauger et al., 
1994) aboard the Hubble Space Telescope. This camera con¬ 
sists of 4 adjacent 800 x 800 Loral CCD’s. Three of these ar¬ 
rays have an image scale of 0.1 arcsec pix“* (WF-chips), while 
the fourth - the planetary camera (PC-chip) - has an image 
scale of 0.046 arcsec pix *, thus well sampling the diffraction 
pattern in the R- and /-band. The FOV in the planetary cam¬ 
era is 25" X 25"and since the jet visible in the ground imagery 
is less than 15 arcsec long, we centered the PC on the posi¬ 
tion of IRS5. We used an orientation that allowed us to detect 
the only visible star within the field-of-view on one of the WF 
chips, leading to potential relative astrometry between epochs 
of ~ 0.1 arcsec. With this orientation we unfortunately could 
not image any of the other brighter nebulosities associated with 
L1551 IR S5. Note that the br ightest object, HH29, has been 
imaged bv IPevine et al.l 1I2OOOI1 . Some fainter nebulosities are, 
however, visible on one of the WF chips. In this paper we only 
discuss features relevant to the jet. 

We observed HH154 during 3 HST visits. During each 
visit, two exposures each were obtained through the F656N 
(Ha) and the F673N (/16717 Ah- 6731A) filters beginning on 
1996.034. Total exposure time for each filter was 2 500 s. On 
1996.078 we obtained 3 exposures each through the F675W (R) 
and F814W (/) filters with the total observing times being 1 800 
s and 2400 s respectively. We again obtained data through the 
F656N (Ha) and F673N (46717 A-h 673 1 A) filters on 1998.074, 
providing us with a time base of almost exactly two years (727 
days). 

The data were reduced and calibrated into units of erg cm~^ 
s ^ us i ng standard methods as d escribed in e.g. Holtzman et al.l 
( I 1995 I) . iHeathcote et akl (Il996l) . iReinurth et al. d 19971) and the 
WFPC2 handbook. 

All the images were registered with each other by using the 
stellar image of VSS4 and the rotation angle information. The 
accuracy in this process is estimated to be ~ 0.1 arcsec. 

2.2. Ground based imaging 

In December 1996 and in October 2001, calibrated images 
were obtained using the Nordic Optical Telescope. The in¬ 
strument combination was identical to that used in the ground 
based spectroscopy except that ALFOSC was set to imaging 


Table 1 . Observed fluxes (erg cm“^ s“') for the assumed work¬ 
ing surface (knot D), the feature F, and the total jet in in all 
of the imaging data. Note that these data have not been de- 
reddened. The quoted accuracy includes both the estimated pre¬ 
cision in the flux calibration and the measurement error. 


Feature 

Ha 

Hp 

[Sii] 


Source 

Accu¬ 
racy (%) 

D 

1.8 x lO-*'* 

_ 

1.8 X 

10-'^ 

FL94 

10 

F 

5.2 X 10-*^ 

- 

3.0 X 

10- 

FL94 

10 

^A //Jet 

3.3 X 10“*'' 

- 

2.9 X 

10-*^ 

FL94 

10 

D 

2.1 X lO-"* 

- 

2.2 X 

10-''* 

HST 1996 

5 

F 

3.8 X 10-’^ 

- 

3.4 X 

10-‘5 

HST 1996 

5 

^A //jet 

5.2 X 10“’^ 

- 

4.8 X 

10-*'* 

HST 1996 

5 

D 

1.8 X 10-'^ 

2.3 X 10“*^ 

- 


NOT 1996 

10 

F 

2.7 X 10-’^ 

2.3 X 10-‘« 

_ 


NOT 1996 

10 

All jet 

4.1 X lO-"* 

5.0 X 10-'^ 

- 


NOT 1996 

10 

D 

1.8 X 10“'^ 

- 

2.0 X 

10-*'* 

HST 1998 

5 

F 

3.4 X 10-’^ 

- 

3.3 X 

10-‘5 

HST 1998 

5 

All jet 

4.3 X 10-’^ 

- 

4.2 X 

10-*'* 

HST 1998 

5 

D 

2.0 X lO-'^' 

- 

1.7 X 

10-*'* 

NTT 1998 

10 

D 

1.8 X lO-"* 

2.4 X 10“*^ 

- 


NOT 2001 

10 

F 

3.8 X 10-'^ 

3.2 X 10-*'^ 

- 


NOT 2001 

10 

^A //Jet 

4.1 X IQ-’'* 

5.2 X 10“*^ 

- 


NOT 2001 

10 


Table 2 . Extinction determined from the Ha/Hy6 ratio, at two 
epochs (December 1996 and October 2001) for several features 
in the jet (see text). The accura cy in the determinatio n of the 
ratios is ~ 20%. Data taken from iFridhinAe7 al]il200-5l) 


Feature 

Ratioj996 

Ratio 2 ooi 

Av,(1996) 

R = 3 

Av(2001) 

R = 3 

D 

8.5 

7.0 

2.8 

2.2 

E 

8.0 

10. 

2.6 

3.2 

F 

11.9 

11.9 

3.7 

3.7 

Region between 

F and end 
of jet (F2) 

12-15 

14-20 

4-4.5 

4.2-5.1 


mode. Images were obtained through a number of filters, of 
which results for Ha, Hy6, [S ii] 4067,68A and [O iii] 5007A are 
reported here, while a more detailed description, as well as the 
presentation of further data is deferred to a later paper. 

2.3. High resolution spectroscopy 

Spectroscopic observations of the jet were carried out with the 
NTT 3.58m telescope - EMMI (ESO Multi-Mode Instrument) 
combination at the European Southern Observatory, La Silla, 
Chile. The observation was carried out on the night of 22 
December 1998. We used the echelle mode with grating No. 
10 and grism No 6, resulting in a wavelength coverage of 
6000A - 8400 A. Using a l"x 10" slit, we have a spec¬ 
tral resolution of 28 000 (equivalent to a velocity resolution of 
~ 11 kms^' at Ha). The seeing was about 1 arcsec, and three 
slit positions, each separated by the same amount in declina- 
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Table 3. Ions detected in spectroscopic observations with the ESO NTT-EMMI instrument. These observations only covered knot 
D (brightest feature in EigureQ. Note that the data are not de-reddened. Note also that the values in the table refer to the maximum 
values. There is a large variation between slit positions for different lines. These differences presumably are due to a geometric 
effect (working surface spatially resolved), discussed in the text. The accuracy in the flux determination is approximately 10%, 
and this error includes both the error in the calibration and the actual measurement error. The error in velocity is 10 km s '. 


Ion 

Rest wavelength 
(A) 

Max line width 
(FWZIkms-') 

Max peak velocity 
(kms^') 

Total flux observed 
(10“'^ erg s“' ) 

Comment 

[01] 

6300.2 

180 

-170 

5.2 


[OI] 

6363.8 

185 

-175 

5.3 


[Nil] 

6548.1 

65 

-185 

0.63 


Hot 

6562.8 

310 

-175 

20.0 


[Nil] 

6583.6 

270 

-185 

7.0 

several separate components' 

[SII] 

6716.4 

200 

-175 

7.26 

£rto(67i7+673i) = 17.2 X 10“'^ erg cm“^ s^' 

[SII] 

6730.8 

220 

-170 

9.9 


[Fell] 

7155.1 

250 

-170 

2.6 


[Call] 

7291.5 

240 

-180 

6.5 

several separate components' 

[Oil] 

7319.4 

- 

- 

- 


[Call] 

7323.9 

125 

-160 

3.0 


[Oil] 

7329.9 

- 

- 

0.16 

Detected at 3cr-level 


' Values in table refer to the strongest component 


tion, were observed. These observations only cover the work¬ 
ing surface of the jet (feature D). Lines detected with the mea¬ 
sured fluxes and velocities are listed in Table |3 The reduction 
of the data were carried out with the ECHELLE context within 
the ESO MIDAS reduction package, and with flux calibration 
being obtained by observing the flux standard HR1544 (see 
http://www.hq.eso.org/observing/standards/spectra for details). 


2.4. Medium resolution spectroscopy 

The second set of spectroscopic observations was obtained 
at the Nordic Optical Telescope (2.6m) in La Palma on the 
nights of 28 & 29 December 1999. We used the ALEOSC im¬ 
ager/spectrograph. This instrument was used in echelle mode 
and was equipped with a Loral 2kx2k CCD with 15 /urn pix¬ 
els corresponding to 0.189 arcsec pixel '. The slit was 13.8 
arcsec long, but due to crowding of the higher orders, only 
12 arcsec has been used. The slit width utilized was 1 arc¬ 
sec. The grating/grism combination has a practical wavelength 
coverage of 4 800A - 11 OOOA, with a spectral resolution of 
4 300 resulting in a velocity resolution of ~ 70 km s ' at Ha. 
Special care was taken in achieving precision pointing of the 
slit on the sky. Thus a short (20 s) exposure through an R- 
filter was obtained at each position. Then without moving 
the telescope, an image of the slit was made on the detec¬ 
tor. Finally, the spectrum was recorded followed by another 
slit image and R-band image before moving to the next pos¬ 
tion. This allows us to spatially register our spectra to within 
one pixel (0.189 arcsec). A total of four (4) positions was 
recorded, with an integration time of 3 600 seconds per po¬ 
sition. The positions of the slits on the sky can be seen in 
Figure |3 Also here the reductions of the data were carried 
out with the ECHELLE context within the ESO MIDAS re¬ 
duction package, and with flux calibration being obtained by 


observing the two flux standards hz4 and BD-i-28°4211 (see 
http://www.hq.eso.org/observing/standards/spectra for details). 

3. Results 

3.1. The wide band HST data 

In Figure ^ we show the R-band representation of the whole 
jet including the location of th e clumps F - D, using d esig- 
nations originally introduced bv iNeckel & Stau^ lll987l) . and 
also used and modified by FL94 and FL98. Then, in Figure|2 
we present part of the /-band images of the jet. The raw images 
have a pixel resolution of 0.046 arcsec per pixel. The HST is 
diffraction limited at the wavelengths in question in this paper, 
and the PC chip is therefore Nyqvist sampling the PSF. The im¬ 
age scale corresponds at the assumed distance of LI551 to ~ 7 
AU or 10'^' cm per pixel. The data are presented at two different 
levels of contrast in order to bring forth both the bright small 
scale structure and the faint extended emission tracing the out¬ 
line of the jet. The R- and /-band data both show a highly struc¬ 
tured jet, but with significant differences. In the R-band we can 
discern several major arc-like components, each consisting of a 
number of knotlike features. Most obvious, a relatively bright, 
narrow and ‘patchy’ structure with a clear ‘wiggle’ ends in a 
very bright feature ~ 10 - 12 arcsec downstream. To conform 
to FL98 we will designate this the northern jet. The second 
component is fainter, has only 1 or 2 knot-like features before 
bending back and intersecting the bright feature at the end of 
the northern jet. Here we designate this feature the southern 
jet. We also see a faint very narrow and short arc (designated 
the northern arc), a few tenths of an arcsec north of the bright 
northern jet (see also Figure|2j- 

The /-band flux in the inner regions of the jet, are enhanced 
relative to the flux in the /?-band image, while the northern jet 
is relatively weaker in the /-band images. All the knots and 
features visible in the /?-band can, however, be detected also in 
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Table 4. Ions detected in medium resolution spectroscopic observations with the Nordic Optical Telescope ALFOSC instrument. 
These observations covered essentially the whole jet (one slit position missing, see Figure 1^. Note also that the values in the 
table refer to the maximum values. There is a large variation between different slit positions for the same ions. These differences 
presumably are caused by the two distinct jets. The total flux refers to the flux from all of the jet (excluding that emanating from 
the missing slit position). Note that the data is not de-reddened. The accuracy in determining the flux is 10%, which includes 
both calibration error and the actual measurement uncertainty. The velocities can be estimated to a precision of ~ 1/2 of a pixel, 
corresponding to about 25 km s'. 


Ion 

^0 

Feature 

FWZI 

Peak Vel. 

F,o, 

Max vel. 

Comment 


(A) 


(kms^') 

(km s“') 

(10“'^ erg s“' ) 

(kms^') 


[Nil] 

6548.1 

D 

200 

-160 

2.5 

-260 




E 

- 

-190 

0.34 

- 




F 

200 

-210 

0.72 

-305 




F2 

300 

- 

0.49 

-400 


Ha 

6562.817 

D 

300 

-150 

17.23 

-310 




E 

270 

-128 

1.43 

-243 




F 

190 

-220 

4.1 

-310 




F2 

160 

-243 

0.82 

-310 




F3 

100 

-50 

0.72 

-105 


[Nil] 

6583.6 

D 

300 

-160 

6.2 

-320 




E 

- 

- 

- 

- 




F 

275 

-250 

2.7 

-345 




F2 

250 

-250 

0.67 

-410 


[SII] 

6716.42 

D 

275 

-215 

5.4 

-325 

S//„.i(6717+6731) = 12.2“ 



E 

- 

- 

- 

- 




F 

180 

-280 

1.35 

-395 

2//i(i(6717+6731) = 3.0" 



F2 

400 

- 

0.5 

-360 

S//m(6717+6731) “ 1-43" 



F3 

180 

-125 

1.49 

-215 

£//i«(6717+6731) = 3.6" 

[SII] 

6730.78 

D 

275 

-160 

6.9 

-325 




E 

180 

-160 

0.54 

-250 




F 

180 

-260 

1.65 

-395 




F2 

400 

-70 

0.93 

-360 




F3 

180 

-115 

2.1 

-215 


[Fell] 

7155.14 

D 

210 

-170 

2.4 

-270 

FWHM= 160 km s '" 



F 

300 

-295 

2.6 

-420 

FWHM= 140 km s ' 



F2 

250 

-330 

0.78 

-430 

FWHM= 150 km s ' 

[Fell] 

7172. 

D 

- 

- 

0.33 

- 


[Call] 

7291.46 

D 

250 

-165 

6.1 

-290 


[Nill] 

7379.6 

D 

- 

- 

1.95 

- 

b 

[Fell] 

7452.54 

D 

- 

- 

0.86 

- 

b 

[Fell] 

8619.33 

D 

250 

-120 

0.64 

-250 

b 


“ erg cm ^ s ' 

''Within a 2" x 2" aperture 


the /-frame. The /?-band filter admits emission from the [Oi] 
6300, 6364A, [Nii] 6548, 6583A, Her and the [Sii] 6717, 

6 731A lines, as well as some scattered light. 

A few emission lines fall within the band-pass of the I- 
hlter viz. [Feii] 7 155A, [Caul 7291A and [On] -H [Can] 

7 325 A - see also lCohen & FullM(ll985h . Although, these lines 
can be seen in our medium and high resolution spectra of the 
jet to be at signiheant flux levels (~ j FjHtr) for feature D) 
for a few of the features, as what concerns the innermost part 
of the jet, they are essentially nonexistent - see Table The 
spatial distribution of this emission - as seen in the spectra - is 
essentially identical to that of the stronger emission lines (e.g. 
[Sii] - see Figure |3l. Features detected in the /-band image 
that have no or little counterpart in the Ha and [S ii] images 
we therefore judge to be primarily due to scattered light. 


3.2. Narrow band HST data 

In Figure|3we display contour plots of our Ha and [S ii] cho¬ 
sen to both bring out faint details and show changes due to e.g. 
proper motions. 

No trace of the ‘northern arc’, so clearly visible in the I - 
frame, can be traced in the emission line images. The ‘southern 
jet’ is visible in the Ha and [S ii] images, but it is found that al¬ 
most all of the line emission originates from the ‘northern jet’. 
Morphologically, the line emission images look fairly similar 
to each other. Major differences first become obvious when we 
analyse the working surface - feature D. In Figure |3 contour 
plots showing the detailed structure of the feature D in the light 
of these two ions, are displayed. The bright core seen in the R- 
and I-band images is obvious in the Ha image, but less so in the 
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[S It] data. Instead in the latter band we see a mottled structure 
with a number of point-like emission centers with the brightest 
emission further away from IRS5 than the bright core seen in 
the Ha data (see also Figure]^. Referring to the ground based 
data of FL94 it was also remarked upon that the [S ii] intensity 
maxima were ‘further out’ fromIRSS. 

Using images obtained under sub-arcsec seeing conditions, 
and then applying an unsharp mask, FL94 found a separation 
of Ha and [S ii] of 0.5 arcsec. This is confirmed by the present 
HST data set. 

3.3. Spectroscopic results 

In Tables 13 and |4] we present a summary of our spectroscopic 
results. The high spectral resolution data obtained towards the 
working surface (feature D), is presented in the hrst table. Here 
the peak radial velocity is found between -160km s ' and - 
185km s'. The maximum FWZI (Full Width at Zero Intensity) 
is ~ 300kms ', which should be ident ihed with the shock ve¬ 
locity at the apex llHartigan et al.Fl 9^ . This is also conhrmed 
by the observation of [O iii] emission from the working surface 
hv lFridhinde7^(l200.5l) . 

In Table |4] we show the medium resolution spectroscopy 
obtained through slits positioned such as to map out the ra¬ 
dial velocity held of the two jets. Also in Figure|3we demon¬ 
strate both the positions where the spectraq were obtained, 
as well as samples of the spectroscopic information (Ha and 
[Fell] 7 155A). These spectra clearly show our conclusions 
drawn about the velocity held in sections 3.4,4.1 and 4.2. From 
the spectroscopic data, it is clear, that there is indeed two ve¬ 
locity systems connected with the jet. A comparison between 
the spectroscopy and the images allow a one-to-one correspon¬ 
dence between the identihed features (F3 - D), and the velocity 
components in the spectra. These results are presented in detail 
in Table|4] The highest radial velocities are found to be associ¬ 
ated with the northern, brighter component, where essentially 
all of the visible shock induced knots are found. A lower (more 
positive) radial velocity system is connected with the southern 
component. The highest velocity found in any ion within the 
jet appear to be associated with the innermost features (F, F2), 
where peak velocities of 430 km s ', and FWZI of up to 400 
km s“' are present (see also section 4.2). A c omparison betwe en 
our data and the IR [Fe ii] spectroscopy of IPvo et al.lll2002i) is 
fully consistent when taking into account differences in sen¬ 
sitivity. A more thoro ugh discussion abou t the spectroscopic 
results can be found in lUiseau et alJj2005al) and iFridlund et al.l 

1 I 2 OO 5 I) . 

3.4. The nature and proper motions of the clumps 
along the jets 

With one major exception, all well dehned features that can be 
identihed in both HST epochs are to be found in the northern 
jet. The exception is the bright oval feature seen in the I-band 
(Figure|2] see also discussion). 

Proper motions of features in the jet were found by 
iNeckel & Staudd (Il987l) . and further studied by FL94. These 


ground based observations led to the identihcation of the fea¬ 
tures denoted A-F. Features A & B are found further out than 
the others and are thus not on the PC chi p of the WFPC2. 
We re fer to FL94 for illustration. In both iNeckel & Staudel 
(Il987h and FL94 it was found that distinct new knots appeared 
from behind the obscuration along the line of sight to IRS5. 
FL94 determined the transverse velocities for knots D, E and F 
to be between 120 km s ' and ~ 300 km s '. Knot E appeared 
well dehned in the middle of the jet in the data from 1989, but 
in 1993 it seemed on the verge of merging with the working 
surface D. At the same time there were indications of a new 
feature appearing at the original spot of E. When comparing 
the separate [S ii] and Ha images, we hnd that the features ap¬ 
pear to be well dehned in each hlterband image, but are found 
at somewhat different positions. This could be explained essen¬ 
tially through a change of excitation conditions at different po¬ 
sitions (see Tableland Eigure|3and|3}- In some of the features 
there are indications of them being representative of clumps of 
gas moving at a high speed t hrough the surroundin g gas (so 
called interstellar bullet (e.g. iNorman ^ '^ll979l) V This is 
especially true for knot E, which is compact and appears at dif¬ 
ferent positions at the two HST epochs. This feature is, how¬ 
ever, currently traveling at transverse velocities about half of 
those measured by EL94. This behavior could be explained by 
the present feature E being either: 

- A different bullet - ejected at a different velocity 

- The feature E observed by EL94 - from the ground - was 

the seeing-smoothed agglomeration of several ’’bullets” 

- A combination of both these options 

As is evident from Table |3 the knots D2 and E2 (see Eigure 
^ are the knots with the highest velocities, consistent with the 
results of EL94. 

Our radial velocity results (see Eigure |3 indicate that the 
highest velocities for a given inclination should be found as 
far in towards IRS5 as possible, which is consistent with E2 
having the highest transverse velocity (see Table |3i- In Eigure 
13 we have zoomed in on the area around knot E. A number 
of smaller, well dehned (sub-)features are obvious in the im¬ 
ages. This F-complex consist of sub-features E2, E’ and E. The 
E2 knot is identihable as a bright isolated knot, of intermedi¬ 
ate excitation, seen in both hlters and in both epochs, while 
E’ has low-excitation characteristics and is harder to identify 
unambiguously in all hlters/epochs. Einally, the E feature is of 
higher excitation. 

The lengthening of the jet with time, as measured by its to¬ 
tal length at each epoch is shown in Table |3 where we have 
collected all the data for this parameter that we have available. 
The average change in length found during this period (1 Dec 
1989 to 21 Oct 2001 — 4 342 days) is 0.3 arcsec yr“' represen¬ 
tative of a transversal velocity of235 kms^^. 

In TableQ we compare the velocity in the plane of the sky 
of the knots A, B, and C (see EL94 for location) over total 
timespan covered by the data with the results in our earlier 
paper. There is maybe an indication that the knots are slow¬ 
ing down as they progress further out beyond the working sur¬ 
face, as well as becoming somewhat less well dehned. This 
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Table 5. Proper motions measured for features in the jet. Note that these data are obtained only from the HST 1996 & 1998 
observations. The actual precision of measurement is 0.1 arcsec. Note that the velocity data for the F-complex is somewhat 
uncertain, since the identification of individual features (particularly F’) have some ambiguity. 


Feature 

Ha 

[Sii] 

Average 

Velocity 

Velocity® 150pc 

Comment 


arcsec 

arcsec 

arcsec 

arcsec yr“‘ 

kms“‘ 


F2 

0'.'76 

0'.'73 

0'.'745 

0''35 

250 

See text about F-complex 

F’ 

0'.'6 

0'.'7 

0'.'65 

0''3 

220 

difficult to identify unambiguously 

F 

0'.'5 

0'.'7 

0'.'6 

0''3 

220 

See text about F-complex 

E 

0'.'53 

0'.'41 

0'.'47 

0'.'22 

160 


D2 

0'.'65 

0'.'62 

0'.'635 

0'.'30 

215 


D1 

0'.'44 

0'.'45 

0'.'445 

0'.'21 

150 


D4 

- 

0''44 

0'.'44 

0''20 

150 


D3 

0'.'36 

- 

0'.'36 

0'.'17 

120 



Table 6. The length of the jet at different epochs. A is the num¬ 
ber of days between successive epochs. The measurement pre¬ 
cision is ~ 0.1 arcsec. 


Filter 

Year 

At 

Days 

Length 

Arcsec 

Rate of change 
Arcsec yr“* 

R 

1989 

0 

10'.'4 

na 

R 

1993 

1364 

11 "0 

0'.'16 

R 

1996 

935 

12'.'24 

0'.'48 

Ha + [S II] 

1998 

727 

12'.'73 

0'.'25 

Ha 

2001 

1316 

14'.'32 

0'.'44 


needs to be conhrmed with observations over a longer time- 
base. Nevertheless, these featur es have now existed as w ell de- 
hned entities since at least 1983 (iNeckel & Staudell987h . In the 
meantime they have moved between 600 and 1 000 AU further 
downstream. The northern jet appears to be wiggly - particu- 

Table 7. The change of distance of features A, B and C 


Feature 

Fridlund & Liseau 119941 

1993-2001 


Arcsec yr“' 

Arcsec yr“* 

A 

0'.'29±0'.'06 

0'.T±0'.T 

B 

0'.'26±0'.'06 

0'.'18±0'.'1 

C 

0'.'26±0'.'05 

0'.'20±0'.T 


larly in the Ha band image. The amplitude of the excursions is 
~ 0.5 arcsec, and the distance between the nodes is 3 arcsec or 
~ 450AU (~ 7x10'^ cm). In unpublished ground based data 
from 1996 (Ha), we have used an unsharp mask on the data, 
and thanks to a much higher signal-to-noise ratio confirmed 
the size and extent of the pattern. 

FL98 interpreted the HST data from 1996 as indicative of at 
least two physically separated velocity components of the jet. 
In our new spectroscopic observations we have complete cov¬ 
erage of the jet with medium and high resolution spectroscopy, 
which allows us to determine the velocity held of the entire jet 
(Figure!^. From different line ratios we can determine some 
of the physical parameters. In order to do so, we use previ¬ 
ously unpublished data, where we have imaged the whole SW 


L1551 complex through narrow interference hlters. Those bi¬ 
ters included Ha, Hy6 and the [S ii] 4067/68A ion. We can then 
a) deredden our data and b) determine the electron tempera¬ 
ture and electron density along the jet, and in and around the 
working surface (knot D). The electron density, n^, can be de¬ 
termined from the [Sii] 6717A/[Sii] 6731A ratio, and a suit- 
able model. We use t he same recent 5-level model for [S ii] as 
iFridlund et al.l(ll998l) . The electron temperature, Tg, can be de¬ 
termined from the same data together with information about 
the hux from the [S ii] 4067/68A. Although inherently faint, 
both knot F and knot D are detected in the blue [S ii] lines, 
and with Huxes of 1.25 x 10“'® erg cm“^ s“' (S/N ~ 4) and 
7.05 X 10“'® erg cm“^ s“' (S/N ~ 11). Note that the apertures 
used are different. Dereddening is achieved from the Ha/Hj6 ra¬ 
tio (see Tables [2 and 13 . Using the above mentioned spectro¬ 
scopic data, and our [S ii] model, we arrive at values of r^(D) 
~ 8 700 K+1 000 K, r,(F) ~ 14250 K+1 750 K, n,(D) ~ 1 600 
cm“^ + 200 cm“^, and ne(F) ~ 1 800 cm“^ + 200 cm“^. Note 
that these values then are averages over the aperture used in 
each case. For knot F this is about 300 x 300 AU^, while for fea¬ 
ture D it is ~530 X 530 AU^, assuming a distance of ~ 150pc. 

We have also analysed our [Fe ii] data using a model atom 
including 142 levels with 1 438 transitions. Our model incor¬ 
porates radiative transfer, but in this case this is unnecessary 
since the optical depth in the [Feii] lines is of order 10“^. We 
have detections of [Feii] 7 155.1 A, 7 172.0A, 7452.5A and 
8 6 I 9 . 3 A for knot D, and our best ht indicates higher and rig 
than the [S ii] data, viz. 13 500K and 3 200 cm“^, respectively. 
Given the somewhat different apertures through which we are 
averaging our data, the values are quite reasonable, however. 
A puzzling result, however, is the absence of the [Fen] line at 
6 875.8A. In our model, which assumes a solar abundance, the 
line at 7454. 6A should always be weaker. We note, however, 
that in observations of HHl, the same situation occ urs, i.e. the 
detec tion of 7454.6A, 8 619.3A, but not 6 875.8A jSolfetalJ 
Il988li . This suggests an explanation. Like in the case of HHl, 
we do not detect any permitted Fell lines. As pointed out by 
ISolf et al.1 lll988l) . the permitted Fell lines require an excita¬ 
tion energy of, on average, 5.2 eV, while the forbidden lines 
have lower excitation energy. All of the detected lines in our 
spectrum have excitation energies of 1.7 - 2.0 eV, while the 
6875. 8A line, falls in between since it requires ~ 3.8 eV. It 
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appears thus that the kinetic energies in knot D are not high 
enough to excite forbidden lines with higher energies than ~ 2.0 
eV. There is no indication of any line with a higher excitation 
energy than that, although a marginal detection of 6875.8A is 
found at the location of knot F, and at a radial velocity of — 
300 kms~^ . An alternative exp lanation is that the atomic data, 
iluinet et alJ lll996l). for the 6875.8A line are in er¬ 


as given in I 
ror. 


3 . 5 . The extinction towards the jet 


In Table [2 we present the flux from some of the persistent fea¬ 
tures found along the jet as a function of time (including some 
of our unpublished data). The presence of HyS emission allows 
us to determine the extinction towards separate parts of the 
jet (assuming the appropriate recombination case), where one 
would expect the extinction to grow as one progresses towards 
the position of IRS5. 

iFridlund et alJ (l2005l) have detected [Out] emission that 
can be taken as direct evidence of a strong shock being present 
at the position of knot D. Tentatively, we then identify the fea¬ 
ture D as a (relatively) strong shock. Somewhat weaker [O iii] 
emission is also detected from the position of knot F. It is at 
these positions we then expect that the collisi onal contribu¬ 
tion t o the population of Ha to be insignificant dHartigan etai] 
Il987ll . We can thus use the Balmer decrement for case B re¬ 
combination in order to estimate the extinction. 


As can be seen from TablesHandlJl in 1996 we find a ratio 
of Ha/HjS of 8.5 for feature D, while the data obtained in 2001 
indicate a ratio of 7.0 for the same object. During the inter¬ 
vening time, this feature has m oved away from I RS5 a t a rate 
of ~ 0.2 arcsec yr ' (Table 1^. ICohen & Fulleil (Il985l) found 
a Ha to Hy6 ratio of ~ 25. This value, however, should be used 
with some caution, since we note that they observed using blind 
offsets from a nearby star, and the Ha and the H/? spectra were 
obtained at different times (pointings). There could thus be an 
offset by an unknown amount between Ha and Hy6. Further, our 
data indicate that the proper motion of this feature over ~ 13 
years, would implicate a displacement of ~ 2 - 3 arcsec, which 
means that the physical conditions could be very different. 


Using a normal galactic extinction curve 
dRieke & Lebofskv 1 19851) our data for knot D indicate 
Ay ~ 3 magnitudes in 1 996 and ~2 magnitudes i n 2001. A 
6 Orionis extinction law dCardelli & Clavto^ll988l) gives val¬ 
ues of 3.6 and 2.9 magnitudes of visual extinction respectively. 
Applying the data of Cohen and Fuller (1985 - data obtained 
1982) we have a visual extinction of 5.8 and 7.4 magnitudes 
depending on the extinction law applied. 


As we progress further in along the jet, towards IRS5, the 
extinction increases. It is close to 4 magnitudes at the position 
of feature F. Towards the very tip of the jets, i.e. the apex where 
they suddenly emerge from invisibility, we find that the ratio 
is indicative of an extinction of 4 to 6 magnitudes. These val¬ 
ues are consistent with the extinction calculated from the X-ray 
spectrum for the L1551IRS5 source dFavata et al.l2002l) . 


4. Discussion 

4.1. The working surface and the physicai conditions 
aiong the jets 

The first, most easily accessed, physical parameter determined 
from our HST data is the level of excitation as d efined by the 
ratio of [S ii] flux to Ha flux jRaga el alJll996ll . For the pur¬ 
pose of this paper we will call this ratio ’’the excitation index”. 
To display the positions where the low excitation [S ii] emis¬ 
sion dominates over the Ha emission - and vice versa we have 
computed the difference images between Ha flux and [S ii] flux 
(e.g. Figurel^. In this figure the most obvious feature is the end 
of the jet (Feature D) where prominent [S ii] , representative 
of low-excitation emission and displayed in white is wrapping 
around the prominent Ha (displayed in black) which is repre¬ 
sentative of high-excitation emission. Since we know the direc¬ 
tion and rough inclination of the jet, the morphology leads to 
the interpretation as being a nice example of a working surface 
with the bow-shock (strong [S ii] ) wrapping itself around the 
Mach-disk (strong Ha). 

Taken together with the excitation index, as defined above, 
the ratio of [O iii] / H/3 refines the understanding of the exci¬ 
tation conditions. According to lRaga et al.l Jl99H) . a value of 
[O III] / Hy6 larger than 0.1 together with a [S ii] /Ha smaller than 
1.5 indicates high excitation, while opp osite values are i ndica- 
tive of very low excitation conditions. IFridlund et alJ (l2005l) 
have found [O iii] 5 007 A emission emanating from both fea¬ 
tures F and D, as well as Hj5 emission distributed in the same 
way as the Ha emission (see discussion about the extinction 
above). 

Feature F is found by IFridlund etalJ ll200.5l) to have an 
[Oiii] /HyS ratio of 0.3+0.06, while the ratio for knot D1 is 
0.22+0.02, both ratios measured in 2001. These are values in¬ 
dicative of a high level of excitation. For the [S ii] /Ha ratio we 
have data covering a larger time span. Feature F (see Figure^ 
has an [Sii] /Ha ratio varying from 0.6+0.1 in 1994 to 1+0.1 
in 1998, thus indicating that an initial high excitation feature is 
somewhat weakening with time. An analysis of the F-complex 
in the HST data shows that the 3 features F2, F and F’ had [S ii] 
/Ha ratios of 1.0+0.1,1.9+0.2 and 0.9+0.1 in 1996. In 1998 the 
values for the ratios were 1+0.1, 0.9+0.1 and 0.8+0.1 respec¬ 
tively. The feature F as observed from the ground, is of course 
the agglomeration of all three subcomponents. These observa¬ 
tions demonstrate the occurrence of changes in the jet not only 
with respect to the morphological properties, but also to the 
physical conditions within a time scale of a few years. Knot F 
is also the second brightest emission feature along the two jets. 
Careful analysis of the HST data show that a very high level 
of excitation is mainly found along the outer (northern) edge 
of feature F, while a very thin (~ a few pixel thick) layer of 
lower excitation (demonstrated by strong [S ii] emission in our 
difference images) is found along the southern edge of the jet. 
Taken together with the [O iii] / HyS index mentioned above we 
conclude that (at least parts of) knot F i s of high but vari able 
excitation in the classification scheme of iRaga et al.l(ll996l) . 

The fainter feature E has an [S ii] /Ha ratio of 0.6 both in 
1996 and in the 1998 data. This is thus indicative of a relatively 
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high level of excitation, regardless of the fact that the feature 
has moved a significant distance i n the intervening time . It is, 
however, not detected in [O iii] bv iFridlund et alJ 1120051) . 

Finally, D, the working surface, is found to have an aver¬ 
age level of excitation index of ~ 1, with very little variation 
over time (see Table If we analyse the separate, very bright 
knots that make up what we have been calling feature D, we 
find a more complex picture. In the Ha and the [S ii] images 
(see Figure 13), we see a number of bright (sub-)features. The 
brightest of these are the ones where we argued in FL98 that 
they were the signature of the working surface of the jet with 
the Mach disk brightest in Ha and the forward shock (or bow- 
shock) being brightest in [S ii] respectively. They are now des¬ 
ignated D1 (brightest in Ha) and D4 (brightest in SII). D1 has 
an excitation index of ~ 0.5 making it a high excitation ob¬ 
ject. This is also supported by the [Out] / HjS ratio of 0.22 (see 
above). The opposite is true for D4 as well as for the ‘bow- 
shape’ surrounding the core. Here we find excitation indeces 
of ~ 1 - 5 indicating low-excitation conditions. 

We already identified the entire feature D as the working 
surface of the jet in FL98. The characteristic features allowing 
us to make this identification are clearly visible in Figure |6l 
Also the interface between the regions dominated by Ha and 
[S It] emission is clearly discernible in Figure 0 which shows 
contour plots of the two HST epochs and where the [S ii] emis¬ 
sion has been subtracted from the Ha emission. The latter emis¬ 
sion is strongly concentrated towards the inner part of the work¬ 
ing surface, at the position of the high excitation core, while 
the [S It] is ‘swept’ into a nice bow-shape surrounding it. The 
steepest gradient of the Ha emission is on the side away from 
IRS5 and on the interface with the region dominated by [S ii] 

. Comparing with the individual images of Ha and [S ii] in 
Figure|3we see that the [S ii] have its steepest gradient where 
the bow-sha pe is furthest aw ay from IRS5. In the context of 
the model o Auartipar^ we interpret these gradients to 

represent the locations of the reverse shock or Mach disk (Ha), 
and the forward, or bow-shock (SII). The bright Ha core of the 
working surface has a FWHM of 0.3 arcsec at right angles to 
the major axis of the outflow, while the whole ‘working sur¬ 
face’, knot D, has a FWHM of ~ 2 arcsec. The separation be¬ 
tween maximum [S ii] (D4) and maximum Ha (Dl) emission 
is 0.6 arcsec or about 2 x 10'^ cm deprojected for a viewing 
angle of 45'’. We are thus resolving the structure of the shock. 

This s epara tion is consistent with the calculations of 
\Hartiear]i \l98^ that — given a cooling distance of ~ Rjet (= 
the jet radius) - show that both the forward shock, i.e. the bow- 
shock, and the Mach disk will be radiative and that the sep¬ 
aration should be just equal to ~ Rjet, as is true in this case. 
In our data both the cooling length - i.e. the distance between 
Dl and D4 - and the jet radius appear to be of the same or¬ 
der (~ lO'^ cm). According to the model of Hartjgan ill989l) . a 
Mach disk which is - as in this case - significantly brighter and 
here also with a much higher le vel of excitation than the bow- 
shock (see lFridlund et al.ll 19971 also Fig. 3) is indicative of a 
less dense jet penetrating into a denser ambient medium. It is 
at this interface where it is possible that the transfer of momen¬ 
tum from the fast atomic wind to the slower molecular outflow 
is taking place. 


The model of lHartig^ lll989l) predicts the relative density 
of the jet and the ambient medium surrounding the jet in terms 
of the relative brightness of the bow-shock and the Mach disk, 
and the jet velocity. For the HH154 jet, we find the density of 
the jet: 

Pjet (0.05 ~ 0.1) pambient (1) 

An analysis taking into account the ng derived in the previ¬ 
ous section and the molecular data of Fridlund & Knee (1993) 
then shows that the ionisation fraction of the jet is ~ 1 (FL98). 
This led to the conclusion in that paper that the optical jet is not 
driving the molecular outflow, since there is a lack of the nec¬ 
essary momentum by at least 2 orders of magnitude. In FL98, 
however, we assumed a velocity of 300 km s ' for the jet mate¬ 
rial (since we then did only have access to low dispersion spec¬ 
troscopy). Inspection of Figure|3 demonstrates that most of the 
material is found at velocities less than this value. Some mate¬ 
rial is traveling at velocities close to 600 km s ', when taking 
the inclination into account, but we note that even if we were to 
allow the bulk of the material to travel at this higher velocity, it 
would still leave us lacking at least one order of magnitude of 
momentum (something which would also require a 100% ef¬ 
ficient mechanism for the transfer of the jet momentum to the 
molecular outflow). In this context, we also note that there are 
no obvious indications of so called ‘turbulent entrainment’, and 
that further the interface for ‘prompt entrainment’ (i.e. at the 
Mach disk) is very small - only ~ 30AU across. In this context, 
we wish to note, that the molecular momentum calculated by 
Fridlund & Knee (1993) refers to the present outflow, i.e. with a 
dynamical time scale of less than ~ a few hundred years which 
is therefore directly comparable to the jet dynamical time scale. 

In the Ha-[S ii] difference images in Figure0we see indi¬ 
cations of a change of position of the Mach disk. If we assume 
that the location of this shock is the sharp interface between the 
[S ii] and the Ha emission, we can see that it has moved (0.5 
arcsec) further out from IRS5 in ~ 2 years. This corresponds 
to a velocity of ~ 180 kms^', assuming a distance of 150pc. 
There is also an indication that the actual apex of the shock is 
turning towards the north - something that could indicate rota¬ 
tion of the jet material. 

4.2. The location of the X-ray source 

Favata et al (2002) discovered 1-2 keV X-rays emanating 
from the jet. Using XMM-Newton, the positional uncertainty 
included most of the jet. Nevertheless, these authors suggested 
that the X-rays were possibly emanating from the working sur¬ 
face in view of it being identified as a strong shock (FL98). 
This was in spite of it having a lower radial velocity and proper 
motion (both ~ 160 kms~') than features found closer to IRS5 
along the iet. lBallv et al.l ll2003l) . using the Chandra X-ray ob¬ 
servatory, have improved on the position, and indeed locate the 
X-ray source in the innermost part of the jet (~1 arcsec from 
the radio position of IRS5). As can be seen from Figure |3 in 
this paper, at this position we find radial velocities for features 
F & F2 of up to 430 km s ', wh ich given the assume d inclina¬ 
tion of the jet of ~ 45'’-60'’ duiseau et al .Il200.5hl) represent 
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true velocities of ~ 500 km s * to ~ 600 km s ^. The immediate 
post-shock temperature is llRagall989l) : 


2.9 X lO^K ^ 

I + X \ 100 km / 


( 2 ) 


Here X is the hydrogen preionisation fraction and an abun¬ 
dance of 0.9 (H) and 0.1 (He) by number has been assumed. As 
mentioned before, FL98 have come to the conclusion, that it is 
very likely that the ionisation fraction is close to unity. 

With these assumptions, we calculate a shock temperature 
of ~ 3 X 10® K to ~ 10^ K , in ex cell ent agreement with the 
X-ray results of iFavata et ai](l2002l) and lBallv et ai1 ( l2003l) . We 
also note here that the material detected at the highest radial ve¬ 
locities probably emanates from the recombination/cooling re¬ 
gion ’downstream’ of the X-ray generating shock, where in all 
likelihood the velocity has dropped somewhat. All this infor¬ 
mation taken together therefore implies that the emission found 
at the highest radial velocities in our data, emanating from in¬ 
side of the feature F2, is associated with the same location (or 
very close to) the source of the X-rays. Thus where the jet first 
becomes visible as it emerges from behind the obscuration. As 
can be seen from Table |2] the extinction toward feature F2 is 
comp atible with the val ues determined from the XMM spec¬ 
trum jFavata et alj2002l) . 

We note here that the presence of an X-ray source at a 
distance of ~ 200 AU from the surface of the disk, as well 
as other shock phenomena generating Lyman photons up to 
maybe 1 000 AU to 1 500 AU further away, are going to have 
an impact on the ionisation structure of the disk surface. This, 
in turn, could have an impact on models trying to explain the 
molecular outflows, such as MHD winds lifting off the disk 
surface. 


4.3. The bright feature in the inner region of the jet 

Both the two jets and the northern arc first appear in the same 
location in the sky (’apex’) and have an opening angle of about 
90°(see Figure |2- the angle between the northern arc and the 
southern jet at their intersection). The jet emission is seen su¬ 
perposed onto a fan-shaped nebulosity that is most evident in 
the I-band filter. The northern arc and the southern jet then bend 
inwards towards the northern jet - the northern arc more dras¬ 
tically than the southern jet - and after 1.4 arcsec all 3 compo¬ 
nents are parallel. Note that it i s within this envel ope that X-ray 
emission is found to originate \Ballv et al. hooA 

Absolute positional information in the HST data is re¬ 
stricted to an accuracy of ~ 1 arcsec, but our images are c onsis- 
tent with the radio position of IRS5 dT.oonev et al .119971) being 
within 0.5 arcsec from the origin of the visible jet. Only 0.6 
arcsec SW from the apex is a very bright feature (see Figure 
|21bottom) visible in the I-band. It is also faintly visible in the 
R-band image, but is actually the brightest feature in the I-band. 
This knot appears elliptical with dimensions of about 0.4 arc¬ 
sec by 0.2 arcsec (FWHM). A careful analysis shows it to con¬ 
sist of two components, with a separation of 0.145 arcsec or 
~ 20AU at the distance of L1551. An approximate R-1 colour 
for both sources is found to be —1-2.0 magnitudes. The fea¬ 
ture is located between the northern jet and the southern jet. 


and at a n angle offset from the outflow direction of those two 
features. Iwhite et ^ll2000l) . have analysed near-infrared HST- 
NICMOS data covering continuum and emission lines between 
\.\2pm and 2.12/im. Inspection of these frames confirm the 
interpretation here, since faint [Fe ii] emission is found from 
knot F, while a strong continuum peak is visible at the same 
position as the oval feature described here. No molecular hy¬ 
drogen emission is detected in the NICMOS data. 

Since the jet emission is mostly caused by recombination- 
presumably through interactions of the outflowing gas through 
shock phenomena, and since the oval feature described here 
is totally absent in the Ha and [S ii] images we suggest that 
it is instead caused by light from the IRS5 system scattered 
up through a polar hole in the circumstellar disk and reflected 
off the inside of the ‘jet-cavity’. Its dimensions of ~ 60AU by 
30AU, with two peaks separated by ~ 20AU, is also consistent 
with the binary nature and par ameters of the IRS5 systern . This 
is supported by the fact that iRodrfgiiez et al.l il2003alhh have 
demonstrated that the jets can be traced in towards the binary 
behind the obscuration, and indeed it thus appears likely that 
they together excavate a ’polar hole’ through which the light 
from the stars/accretion disks can escape. 


5. Conclusions 

In summary, we And; 

- We suggest that the bright, oval feature detected close to 
the apex of the jet in the I-band is caused by light from the 
IRS5 binary system scattered up through a (polar?) hole 
in the circumstellar disk and reflected off the inside of the 
‘jet-cavity’towards our direction. This is supported by the 
colour of the feature. 

- We And that the position of the origin of the X-rays em¬ 
anating from the jet is coincident or very nearly so with 
the location of the observed maximum radial velocities ve¬ 
locities (representative of true velocities of ~ 500 kms^'- 
600 kms^'). 

- We have mapped out the proper motions resulting in trans¬ 
verse velocities, as well as mapping the complete radial ve¬ 
locity field of the jet. A few jet features (e.g. knot E) now 
have very much different velocities. We suggest this is be¬ 
cause the old feature has disappeared and we are actually 
observing a new ’’interstellar bullet” traveling down the jet 
towards the working surface. The time scale for the ejection 
of bullet like features is thus of the order of a few years. 

- The average rate of change of the length of the whole jet, 
found during our observing period (1 Dec 1989 to 21 Oct 
2001 = 4342 days) is 0.3 arcsec representative of a 
transverse velocity of 235 km s“'. 

- The time scale for variations of physical conditions within 
or near the jet (like the excitation conditions or the extinc¬ 
tion) is of the order of a few years. This is similar to the 
time scale for morphological changes. 

- Within the F-complex of knots, we find a number of 
well defined moving features (F, F’ and F2), with large 
proper motions, but with much different excitation condi¬ 
tions among them. 
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Further, we also find; 

- The excitation conditions and the morphology within the 
feature D identify it as a ’’working surface” where the jet 
impacts either the ambient medium or material moving at a 
slower speed (from a previous ejection event). 

- In the latter context, it is noteworthy, that feature D has 
maintained its individual character during more than 20 
years and while traveling more than 1000 AU further 
downstream. 

- The Mach disk and the forward (or bow-)shocks are well 
resolved. The separation between high excitation (as traced 
by Ha), and low-excitation (as traced by [S ii] in the work¬ 
ing surface (knot D) is ~ 0.6 arcsec. 

- The Mach disk is moving ’downstream’ with a transverse 
velocity of ~ 180 km s'. During the time interval between 
our HST observations it also appears to have changed di¬ 
rection. 

- The cooling distance in the working surface of the northern 
jet is ~ equal to the jet radius, Rjet ~ 0.6 arcsec 

Acknowledgements. We gratefully acknowledge valuable discussions 
with Drs. G. Olofsson, A. Heras, T. Prusti and G. Pilbratt. Some of the 
data presented here have been taken using ALFOSC, which is owned 
by the Instituto de Astrofisica de Andalucia (lAA) and operated at 
the Nordic Optical Telescope under agreement between lAA and the 
NBIfAFG of the Astronomical Observatory of Copenhagen. The com¬ 
ments and positve criticism of the referee. Dr. Bo Reipurth, which led 
to a more stringent paper, are gratefully acknowledged. 

References 

Bally, J., Feigelson, E., & Reipurth, B. 2003, ApJ, 584, 843 
Bieging, J. H. & Cohen, M. 1985, ApJ, 289, L5 
Campbell, B., Persson, S. E., Strom, S. E., & Grasdalen, G. L. 
1988, AJ, 95, 1173 

Cardelli, J. A. & Clayton, G. C. 1988, AJ, 95, 516 
Cohen, M., Bieging, J. H., & Schwartz, P. R. 1982, ApJ, 253, 
707 

Cohen, M. & Puller, G. A. 1985, ApJ, 296, 620 
Cudworth, K. M. & Herbig, G. 1979, AJ, 84, 548 
Devine, D., Bally, J., Reipurth, B., Stocke, J., & Morse, J. 2000, 
ApJ, 540, L57 

Devine, D., Reipurth, B., & Bally, J. 1999, AJ, 118, 972 
Pavata, E, Pridlund, C. V. M., Micela, G., Sciortino, S., & Kaas, 
A. A. 2002, A&A, 386, 204 

Pridlund, C. V. M., Bergman, R, White, G. J., Pilbratt, G. L., & 
Tauber, J. A. 2002, A&A, 382, 573 
Pridlund, C. V. M., Djupvik, A. A., Liseau, R., et al. 2005, in 
preparation 

Pridlund, C. V. M. & Piseau, R. 1994, A&A, 292, 631 
Pridlund, C. V. M. & Piseau, R. 1998, ApJ, 499, P75 
Pridlund, C. V. M., Piseau, R., & Gullbring, E. 1998, A&A, 
330, 327 

Pridlund, C. V. M., Piseau, R., & Perryman, M. A. C. 1993, 
A&A, 273, 601 

Pridlund, C. V. M., Nordh, H. P., van Duinen, R. J., Aalders, J. 

W. G., & Sargent, A. I. 1980, A&A, 91, PI 
Pridlund, C. V. M., Sandqvist, A., Nordh, H. P., & Olofsson, 
G. 1984, A&A, 137, P17 


Pridlund, M., Huldtgren, M., & Piseau, R. 1997, in lAU Symp. 
182: Herbig-Haro Plows and the Birth of Stars, Vol. 182, 
19-28 

Hartigan, P. 1989, ApJ, 339, 987 

Hartigan, R, Raymond, J., & Hartmann, P. 1987, ApJ, 316, 323 
Heathcote, S., Morse, J. A., Hartigan, R, et al. 1996, AJ, 112, 
1141 

Holtzman, J. A., Burrows, C. J., Casertano, S., et al. 1995, 
PASP, 107, 1065 

Keene, J. & Masson, C. R. 1990, ApJ, 355, 635 
Kenyon, S. J., Dobrzycka, D., & Hartmann, P. 1994, AJ, 108, 
1872 

Padd, E. P, Puller, G. A., Padman, R., Myers, P. C., & Adams, 
PC. 1995, ApJ, 439, 771 

Piseau, R., Pridlund, C. V. M., & Parsson, B. 2005a, ApJ, 619, 
959 

Piseau, R., Pridlund, C. V. M., & Parsson, B. 2005b, 
Astrophysical Journal in press 

Pooney, P. W., Mundy, P. G., & Welch, W. J. 1997, ApJ, 484, 
P157 

Mundt, R. & Pried, J. W. 1983, ApJ, 274, P83 
Mundt, R., Ray, T. R, & Raga, A. C. 1991, A&A, 252, 740 -h 
N eckel, T. & Staude, H. J. 1987, ApJ, 322, P27 
Norman, C. & Silk, J. 1979, ApJ, 228, 197 
Pyo, T., Hayashi, M., Kobayashi, N., et al. 2002, ApJ, 570, 724 
Pyo, T., Hayashi, M., Kobayashi, N., et al. 2005, ApJ, 618, 817 
Quinet, P, Pe Dourneuf, M., & Zeippen, C. J. 1996, A&AS, 
120, 361 

Raga, A. C. 1989, in Pow Mass Star Pormation and Pre- 
main Sequence Objects, ed. B. Reipurth (European Southern 
Observatory, Garching bei Munchen), 281 
Raga, A. C., Bohm, K. H., & Canto, J. 1996, Revista Mexicana 
de Astronomia y Astrofisica, 32, 161 
Raga, A. C., Riera, A., Masciadri, E., et al. 2004, AJ, 127,1081 
Reipurth, B. 1999, A general catalogue of Herbig-Haro objects, 
electronically published via World Wide Web at http; 
casa.colorado.edu/hhcat 

Reipurth, B., Hartigan, R, Heathcote, S., Morse, J. A., & Bally, 
J. 1997, AJ, 114, 757 

Rieke, G. H. & Pebofsky, M. J. 1985, ApJ, 288, 618 
Rodriguez, P. P, Canto, J., Torrefies, J. M., & Ho, P. T. P. 1986, 
ApJ, 301, P25 

Rodriguez, P. P, Curiel, S., Canto, J., et al. 2003a, ApJ, 583, 
330 

Rodriguez, P. R, Porras, A., Claussen, M. J., et al. 2003b, ApJ, 
586,P137 

Sarcander, M., Neckel, T., & Elsaesser, H. 1985, ApJ, 288, P51 
Snell, R. P., Bally, J., Strom, S. E., & Strom, K. M. 1985, ApJ, 
290, 587 

Snell, R. P., Poren, R. B., & Plambeck, R. P. 1980, ApJ, 239, 
P17 

Solf, J., Bohm, K. H., & Raga, A. 1988, ApJ, 334, 229 
White, G. J., Piseau, R., Men’shchikov, A. B., et al. 2000, 
A&A, 364, 741 






C.V.M. Fridlund et al.: HST and spectroscopic observations of the L1551 IRS5 jets (HH154) 


13 


northern arc 



bright feature 


northern jet 


southern jet 



Fig. 2. The inner part of the 1996 HST I-band image displayed as both a contour plot and a grey scale in order to show both faint 
structures and the details in the ’’bright feature (see text). The 0",0"position is the nominal L1551 IRS5 position. In the bottom 
panel, the lowest contour value is 5.0 x 10“'^ erg cm“^ s“' , and each contour is separated by 1.4 x 10“'^ erg cm“^ s“* . This data 
has been smoothed by a Gaussian filter with a 2 pixel radius. 
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Fig. 3. The 1996 (left upper panel) and 1998 (right upper panel) Ha images displayed as contour plots. The lowest contour is at 
the 3cr level, and each consequtive contour is at intervals of Icr. The 0",0"position is the same as in Figures 1 and 2. Note that the 
innermost part of the jet (position F3 in the R-band image) is not seen in emission. The bright spike in the top panel (7"3,-4'.'8 is 
a Cosmic ray hit.) The 1996 [S ii] image (bottom left panel) and 1998 [S ii] (bottom right panel) imagesare also displayed as 
contour plots with the same contours. The bright pointlike knot in the lower left panel at position 11 '.'5,-2'.'5 is a cosmic ray strike 
that could not be removed. All images are obtained with the same exposure time. 
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Fig. 5. The 1996 and 1998 Ha (top) and [S ii] (bottom) contour representations of feature D - the presumed working surface. 
North is up and East towards the left. The lowest contour is 2.5x 10^'^erg cm“^ s“^ pixeE', and each contour is separated by 
l.Ox 10 '*erg cm“^ s ' pixel '. The bright, compact feature, dominating the Ha frame has been given the designation D1 in 
this paper. Note the proper motions over the 770 days between the 2 epochs representative of transverse velocities of about 150 
kms *. CR is a cosmic ray strike. 
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Fig. 6. An Ha -[S ii] 6717+6731A difference image (1996 data) of the complete jet, shown as a greyscale. Here Ha is represented 
by black, while the [S ii] is white. The picture is smoothed to 0.1 arcsec resolution, and the image is 11 arcsec by 6 arcsec. The 
strong white/black feature is knot D - the working surface. The spatial separation of black (Ha) and white([Sii] ) is interpreted 
in the text as the resolving of the shock into the Mach disk (black) and forward- or bow-shock (white). Note that apart from the 
working surface the rest of the jet cancels out almost perfectly, leaving only a few filaments of separated emission. 



Fig. 7. Difference images of the working surface (feature D - see text) between the Ha and [S ii] images in 1996 (top) and 1998 
(bottom). The positions are relative to an arbitary position. 
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Fig. 8. The Ha and [Fe ii] 7155 A spectra obtained with ALFOSC. In the top row, the different slit positions are shown superposed 
on 20s R-band images. The effective slit width is 12" long. In the middle row are the [Feii] spectra and in the bottom rows are 
the Ha spectra. All are displayed as contour maps, where the contours have been selected to show the maximum extent of 
emission. A radial velocity scale w.r.t. barycentric zero has been used. Comparing on the one hand the spectra of Ha and [Fe ii] 
7155A spectra and on the other the slit positions with repsect to the jet, clearly demonstrate the two different velocity systems, 
originally discovered by FL98. At position 12.2", in both spectral lines, feature D is prominent. A peak velocity of 160 km s 'is 
indicated. Features D, F and F2 have been marked according to their 1998 HST positions 









